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We reduce complex stripped patterns to a basic topological network of edges and vertices to 
define defects and measure their influence on the pattern. We present statistics on the spatial 
and temporal distribution of defects within the state of spiral defect chaos state in experiments 
on Rayleigh Benard convection. These measure the role of boundary influence on the dynamics, 
and suggest an exponential distribution for the length of edges in the pattern. We also indicate a 
systematic method to study hierarchies of defect interactions based on the network structure. 

PACS numbers: 



Topological variations represent basic changes in the 
connectivity within the geometric structure of a system. 
The mechanisms leading to topological change are of- 
ten associated with physically significant events under- 
lying the dynamics. The reduction of patterns to their 
discrete topological features is therefore useful in build- 
ing low dimensional phenomenological models for the dy- 
namics. In some systems, a transition between different 
geometric structures occurs at spatially and temporally 
localized regions. For example, the patterns exhibited 
in Rayleigh-Benard convection are observed to transition 
from regular structures of parallel rolls or target patterns 
to complex time-dependent chaos through the formation 
of localized defects These defects mark regions of 
merger and pinch-off of the stripped patterns exhibited 
resulting in changes in the topology Understanding 
the dynamics of these defects is therefore important to 
understanding the irregular and chaotic patterns formed. 

Defects highlight regions where symmetries of a pat- 
tern are broken. The identification of a defect therefore 
often relies on measures of symmetry. In regular convec- 
tive patterns, defects have been identified by marking the 
variation in the phase of an underlying parallel striped 
pattern 0, H, . Spatio-temporal chaotic systems often 
exhibit intricate patterns with interwoven structure. The 
absence of symmetries in these patterns makes it difficult 
to define and identify defects within the structure. In this 
paper, we present a technique to reduce many conven- 
tionally observed complex striped patterns to a discrete 
set of vertices and edges representing a graph network. 
This allows fluid systems to be studied using many of the 
techniques developed in studying complex networks. 

The patterns of fluid flow during Rayleigh-Benard con- 
vection represent among the most thoroughly studied 
patterns in physics |6|, |7|. In this system, a thin layer 
of fluid (with depth d) is held in a vertical thermal gradi- 
ent (with temperature difference AT). At a critical tem- 
perature difference, ATc, buoyant forces overcome ther- 
mal and viscous dissipation to result in a convective in- 
stability. The temperature difference across the cell is 
parameterized by the non-dimensional reduced Rayleigh 
number, e = (AT - ATc)/ ATc). Ideally, the velocity of 
the fluid is in the form of a set of parallel rolls. A two 
dimensional projection of this structure along the hori- 




FIG. 1: The above sequence of images illustrates the two 
main steps of the image analysis. In (a) a unfiltered image 
of the convection rolls during spiral defect chaos is shown. 
This image is filtered and thresholded to yield a binary image 
with two components of equal area corresponding to bright 
and dark regions as seen in (b). Bright and dark regions 
are independently reduced to one dimensional curves, while 
maintaining the topology of the original pattern. The final 
skeletonized image is shown in (c). The vertices in the bright 
and dark regions have been marked by circles. 



zontal direction yields a parallel stripped pattern. The 
general, convection rolls form stripes that may be sim- 
ply or multiply connected regions. A transition from a 
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stationary and regular parallel striped pattern to a time 
dependent chaotic state known as spiral defect chaos has 
been observed with an increase in thermal driving 

Our experiments use gaseous CO2 pressurized at 
30atm as the convective fluid confined within a circular 
convection cell of diameter 3.8cm and d = 0.069cm. We 
study the patterns exhibited during spiral defect chaos 
at e = 1.1 to demonstrate the techniques presented in 
this paper. This state is known as a classic example of 
spatio-temporal chaos. The temperature difference be- 
tween hot upfiows and cooler downfiows are also associ- 
ated with a variation in refractive index of the gas. These 
patterns are observed as bright and dark regions in shad- 
owgraphic images [9]. In the experiments, the images are 
captured using a digital CCD camera with a resolution of 
650 pixels x 515 pixels. Each gray-scale image obtained 
is Fourier filtered to remove optical distortions inherent 
in the measurement process (pixel noise within the CCD 
array for example) and cropped to the circular region of 
the convection cell. The filtered images are thresholded 
to their median value of intensity to produce a binary 
image. This threshold reduces the image to two com- 
ponents with equal areas representing bright and dark 
regions in the original image (see figure [TJd). 

We process each component of the binary image sep- 
arately. This is done by eroding the pixels at the 
boundary, except when such erosion would result in con- 
nected regions breaking apart. Repeated iteration of the 
procedure on the components of the binary images re- 
duces their structure to be one pixel thick. Some spu- 
rious structures also result with the formation of nu- 
merous spurs. These spurs result from a variation in 
the thickness/ wavelength of the convective stripped pat- 
terns. These artifacts are removed by further eroding 
terminating pixels of the structure. This erosion is done 
to correspond to approximately half the wavelength of 
the stripped pattern. The final result of the procedure 
is illustrated in figure [T] where the raw experimental im- 
age(a), the thresholded image(b) as well as the skele- 
tonized image(c) are indicated. Further details about 
this morphological procedure is found in [10]. 

A limitation inherent to the technique described is that 
the spatial resolution of the resulting skeletonized fea- 
tures is one pixel. This limitation arises because the 
skeletonization technique applies to binary images. Sub- 
pixel spatial resolution may be obtained by weighting 
the localized regions around desired features by intensity 
variations in the original gray-scale images. The results 
reported here are robust within the current limits on spa- 
tial resolution. 

The image analysis procedure outlined provides an al- 
gorithm to reduce two dimensional stripes to one di- 
mensional lines. This skeletonized description maintains 
the basic topological structure of the pattern by preserv- 
ing the Euler number. By construction, the one dimen- 
sional curves representing bright and dark regions do not 
intersect each other, however, they often self-intersect. 
The intersections form vertices along the one dimensional 



curves and are marked by circles in Figure [TJ^. These 
vertices represent regions of defect formation in the pat- 
tern. The statistics of defects may therefore be measured 
through those of the vertices. The shadowgraph images 
are thus rendered into a network of vertices (defects) and 
edges (convective rolls) between two defects. Lines with- 
out defects also occur in the form of loops or with termi- 
nation. For the purposes of the current study, we refer 
to all line segments as edges. 

The occurrence of only three fold vertices within the 
pattern places a severe constraint on the structure of the 
network. Similar constraints are known in other two- 
dimensional systems such as foam where local surface 
tension forces dictate the equilibrium structures or in reg- 
ular hexagonal packing. In the present case, this tri-fold 
vertex structure arises from the hydrodynamic stability 
of the structure [ll| . 
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FIG. 2: The above scatter plot indicates the positions of the 
defects at different times for bright (a) and dark(b) convection 
rolls. The colormap indicates the time at which the defects 
were detected scaled by the vertical diffusion time. 

During the flow, the spatial locations of the defects 
change with time. Figure [2] indicates the variation in 
position of defects observed within the convection cell. 
Often, the defects are seen to glide along the convection 
rolls. The streaks observed in the scatter plot in Figure [2] 
are a result of such motion. We also observe that the 
number of defects is higher for dark regions. This may 
be because the flow is non-Boussinesq at the parame- 
ter values used. The non-Boussinesq parameter (deflned 
in [12]) has a value of approximately 0.54 close to the 
onset of convection indicating a break in the Boussinesq 
symmetry. While the number of defects shows a break in 
the Boussinesq symmetry, the probability distribution of 
the defects is symmetric as indicated in Figure [3l 

It is diflicult to measure the influence of the boundaries 
in driving the formation of defects. The radial distribu- 
tion of defect positions is shown in Figure [3l The square 
of the radius is used along the x-axis as the circular area 
considered varies quadratically with the radius. The flg- 
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FIG. 3: The radial distribution of the defect density is shown 
above. The x-axis is normahzed by the radius of the convec- 
tion cell. The solid and hollow markers represent bright and 
dark regions of the image. 



ure shows that the density of defects increases away from 
the boundaries. Further we calculate that the mean ra- 
dial velocity of the defects vanishes independent of their 
radial position, suggesting that the flux of defects moving 
radially inward is balanced by those moving radially out- 
ward. The probability distribution in Figure [3] is there- 
fore stationary. The stationarity of the distribution im- 
plies there is no radially localized source/sink for defect 
creation/annihilation. We also note that at the aspect 
ratio studied, the radial distribution of defects does not 
attain a spatially uniform value. It would be interesting 
to study the aspect ratio at which a region with spatially 
uniform defect density is obtained. This region would 
reflect the defect dynamics inherent to the flow during 
spiral defect chaos independent of the boundaries. 

In the state of spiral defect chaos, the motion of the 
convection rolls result in defects being continually created 
and annihilated. The vertices in the network constructed 
are created either through the formation of a new edge 
along an existing edge, or by the intersection of two sepa- 
rate edges. A vertex may be annihilated through the con- 
traction of one of the three edges deflning it. When two 
edges are shared between two vertices, the contraction 
of the shared edges results in the annihilation of a pair 
of vertices/defects. Similarly, progressively higher order 
interactions may be deflned, for example, with a vertex 
triplet, with each vertex sharing a pair of edges with to 
other vertices(as in a triangle) and so on. The different 
processes of edge creation and annihilation occur over a 
range of timescales. We use a standard Particle Image 
Velocimetry program to identify and track the motion of 
individual defects during the flow [l^]. 

The lifetime of a defect varies inversely with the num- 
ber of defects as indicated in flgure [H The temporal 
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FIG. 4: The above plot indicates the amount of time each sin- 
gle defect was tracked in units of the vertical thermal diffusion 
time. The horizontal axis is scaled by the vertical diffusion 
time. The solid lines are fits to logarithm of the probability us- 
ing third order polynomials. The coefficients for 
are (0.002,0.1,-1.48,-1.68) for the solid markers (bright regions) 
and (-0.014,0.39,-3.38,-1.46) for the hollow markers(dark re- 
gions) respectively. 



resolution of image capture in the experiments sets the 
lower bound of the lifetime measurements. We flnd that 
the logarithmic distribution of lifetimes is well flt by a 
cubic polynomial. The coefficient for the linear term for 
this fit is large in comparison to those of the quadratic 
and cubic term which represent small corrections. Also, 
the number of defects with the same lifetime is larger for 
bright in comparison to dark stripes in the images. 

The formation of defects in Rayleigh-Benard convec- 
tion has been associated with instability mechanisms 
about the straight roll state. In our system, the skew- 
varicose instability and cross-roll instability primarily 
drive variations in the pattern at short and long wave- 
lengths respectively [6]. The vertex description of a de- 
fect does not directly provide information on these in- 
stability mechanisms, however, it might be possible to 
observe differences in the network structure around the 
neighborhood of defect formation to suggest the type of 
instability mechanism involved locally. 

Each edge in the network is either interspersed with 
vertices, terminates at the boundary/bulk or very rarely 
forms a loop with no vertices. The edges in the graph 
network represent these various states of the convection 
rolls. The distribution of lengths of the edges is shown in 
figure [5l We find that the number of edges of a given 
length show a dominantly exponential decay with in- 
creasing length of the edges. The deviation from this 
behavior for edges with large lengths may be due to the fi- 
nite diameter of the convection cell. The curvature of the 
convection rolls allow for edges of length exceeding the 
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FIG. 5: The distribution of the edge lengths in the net- 
work structure is displayed in the above figure for bright (soUd 
markers) and dark regions (hollow markers) of the convection 
pattern. A third order polynomial fit (solid line) to the data 
yields the coefficients to the cubic, quadratic, linear and con- 
stant term as (-0.1,2.18,-8.75,-1.64) respectively. The domi- 
nant linear term suggests that the edge lengths follow a pri- 
marily exponential distribution. The lengths of the edges have 
been scaled by the diameter of the convection cell. 

diameter of the cell. Since the width of bright and dark 
stripped patterns are similar, and their areas are equal, 
this also implies that the total length of the edges are 
equal. Figure [5] indicates that the sectioning of lengths 
through defect formation or termination for bright and 
dark regions are equivalent. 



There have been many extensive investigations of de- 
fects, their identification and their role in fluid flow. 
Some of the interesting studies and advances have been 
cited within the literature referred to in this paper. Of 
particular interest are recent studies that report mea- 
sures similar to those extracted in this paper from defect 
trajectories [14]. Geometric diagnostics of chaotic pat- 
terns resulting from defect dynamics too have been re- 
ported [15]. In addition, the study of patterns through 
different topological characterizations has been of inter- 
est. Extracting rigorous and dimension independent fea- 
tures from pixelated images however can be daunting. 
The use of computational presents one of the simplest 
and most easily accessible topological characterization of 
patterns [1^. 

Our present study describes a technique to obtain a 
graph network for images corresponding to experimen- 
tally obtained patterns of chaotic fluid flows. These pat- 
terns are similar to those obtained from a number of other 
canonical systems exhibiting spatio-temporal chaos such 
as reaction-diffusion systems, granular systems, electro- 
convection in liquid crystals etc. (See Q for a review). 
Applying techniques described here lay a common for- 
malism to describe these systems well into the chaotic 
regime. In the study presented, the patterns are reduced 
to graphs representing two components of the flow. The 
vertices of the graph are identified as defects within the 
pattern. A number of statistics describing the distribu- 
tion of defects and their influence on the underlying net- 
work structure have been described. We have measured 
the lifetime of a defect and the dispersion of convection 
roll lengths. We have also quantified the influence of the 
boundaries on the spatial distribution of the defect dy- 
namics. 
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